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ABSTRACT 



We present a re-analysis of our H- and .fT-band photometry and light-curves 
for GCIRS 16SW, a regular peri odic source near th e Galactic center. These 
data include those presented by IPePoy et al.l (120041 ) : we correct a sign error 
in their reduction, finding GCIRS 16SW to be an eclipsing binary with no color 
variations. We find the system to be an equal mass overcontact binary (both stars 
overfilling their Roche lobes) in a circular orbit with a period P = 19.4513 days, 
an inclination angle i = 71°. This confirms and strengthens the findings of 
Martins et al.l (120061 ) that GCIRS 16SW is an eclipsing binary composed of two 
~ 50M© stars, further supporting evidence of recent star formation very close to 
the Galactic center. Finally, the calculated luminosity of each component is close 
to the Edding t on lum inosity, implying that the temperature of 24400 K given by 
Najarro et al.l (119971 ) might be overestimated for these evolved stars. 



Subject headings: Galaxy: center 
binaries: eclipsing 



stars: individual (GCIRS 16SW) — stars: 



Introduction 



Standard star formation modes are thought to break down near a supermassive black 
hole (SMBH), raising the question o f whether or not star formati on near a SMBH is possible, 
and if so, through what mechanism (INayakshin fc Sunyaevll2005l ). Our own Galaxy provides 
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us with a unique opportunity to study individual stars in the presence of a SMBH, namely, 
Sgr A*. Direct observations of massive, and therefore young, sta rs close to Sgr A* ind icate 
that there has been recent star formation at the Galactic center (ILebofsky et al.lll982l ). 



GCIRS 16SW (hereafter IRS16SW) is a variable sou rce near the Galactic center (a = 
17 h 45 m 40 s lj g = _29°00'29", J2000.0). btt et al.l Jl999h reported that the source is regu- 
larly variable and suggested that it could be a binary star with very massive components. 
DePoy et al.l (12004 ) confirmed the period of the object, but argued t hat the source was more 
likely a pulsating variable. Recently, however, iMartins et al.l (120061 ) reported spectroscopic 
observations of IRS16SW that showed radial velocity variations consistent with a binary 
composed of two massive stars. 



Prompted by the convincing nature of the radi al velocity var i ations seen by lMartins et al 



(120061 ) we have re-analyzed the data presented by iDePoy et al.l (120041 ) as well as additional 
data from the same observing campaign. We find that the original data reduction process 
was seriously flawed. In partic ular, the colo r variation, light-curve asymmetry, and sign 
of the brightness variations that IDePoy et al.l presented are artifacts of the data reduction 
process. 



In this Letter, we report on the re-analysis of the IDePoy et al.l data. We find that 
there is no color change in IRS16SW over its variations and that the shape of the light 
curve is consi s tent w ith an eclipsing binary system. The new results are consistent with 
Martins et al.l (120061 ) and confirm that IRS16SW is a binary composed of two massive stars. 
In £j2]we describe the observations and present the data, in §3] we describe the best-fit model 
to the light-curve, and in §4] we discuss and summarize the results. 



Observations 



Observations of the Galactic center in the H (1.6/xm) and K (2.2/im) bands were made at 
the Cerro-Tololo Inter- American Observa tory (CTIO)/Yale 1-m telescope using the facility 
optical/infrared imager (ANDICAM; see bePov et alil2003l for details). ANDICAM has a 
pixel scale of / /22 pix -1 on a 1024 x 1024 array. Both H- and .fT-band images were taken 
in the 2001 and 2 002 observing seasons; if-band data were also obtained in 2000. (The 
DePoy et al.ll2004l analysis includes only the 2001 data.) The observing campaign consists of 
every usable night from UTC 2000 August 13 (HJD 2451769.5) through UTC 2000 October 
14 (HJD 2451831.5), UTC 2001 May 20 (HJD 2452049.5) through UTC 2001 November 3 
(HJD 2452216.5), and UTC 2002 June 9 (HJD 2452434.5) through UTC 2002 September 25 
(HJD 2452542.5). Each night, a set of seven slightly offset images were obtained and then 
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combined and trimmed to form a final nightly image. The if-band images consist of 30 s 
exposures, and it took about four minutes to obtain the group of seven images; the iT-band 
images consist of 10 s exposures and took about two minutes to obtain. 

The final 512 x 512 pixel images, corresponding to a field of view of 112 x 112 arcseconds, 
are approximately centered on the Galactic center. After image quality cuts were made, there 
are a total of 144 if-band and 137 iT-band images. The seeing ranges from 0'/93 to l'/93 
full- width half maximum; in general, the if -band images are of higher quality (typical seeing 
~ l'/3) than the i^-band (with typical seeing ~ l'/45). 

Bec ause the field is crowded, we reduc ed the data using the ISIS difference image analysis 
package (lAlardll2000l ; lHartman et al.ll2004l ) . This analysis revealed a sign error in the original 
DePoy et al.l (120041 ) reduction; IRS16SW is clearly an eclipsing binary. Because IRS16SW 
i s sub ject to significant blending — there are roughly half a dozen sources in the lOtt et al. 



(119991) catalog withi n about one arcsecond of IRS16SW — -we cali brated our light-curv e s with 



the lOtt et al.l (119991 ) data as presented by iMartins et al.l (120061 ). The IMartins et al. 
.fT-band photometry gives a mean magnitude of 0.2 mag higher than the lOtt et al. 
mean, and includes two more seasons of data. 



<S) 

Jl999h 



Ott et al.l found a variability amplitude of 0.55 mag; using DAOPhot pho tometry to scale 



the I SIS fluxes, we find an amplitude of ~ 0.35 mag in both H and K (see lHartman et al. 



2004 . Appendix B). This substantial amplitude difference is indic ative of significant b lend- 
ing in our da t a. W e used the period of 19.45 days reported by iDePoy et al.l (120041 ) and 
Martins et al.l (I2006T) to scale ou r iT-band light-curve to have the same mean magnitude and 
amplitude as the IMartins et aLl data. There were 110 nights for which data were obtained in 
both the H- and iT-bands, providing contemporaneous measurements of the H — K color. 
Using the DAOPhot photometry, we find a constant H — K color with an rms of 0.05 mag; 
this color does not vary with time, phase, H-, or i^-band magnitude, as shown in Figure [TJ 
Lacking properly calibrated if-band data, we scaled the similarly blended ii-band data to 
have the same amplitude and mean magnitude as the iT-band data. These scaled light-curves 
form the basis of our analysis; Table [1] gives the final scaled time-series photometry. 



3. Light Curve Analysis 



We simultaneously fit the 144 if -band and 184 fT-band points (including the IMartins et al. 



fT-band points, except for two noisy points at HJD 2498704 and 2499908) using the October 



2005 version of the Wilson-Devinney (WD) code (IWilson fc Devinneyl Il97ll ; IWilsonl 11979 
1990h in the overcontact mode (MODE 3). We fixed T eS1 = 24400 K as estimated by 
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Najarro et al.l (119971 ) and used th e squar e -root limb darkening law, ta king the value s of the 
limb darkening coefficients from I Clare tl (120001 ) for a LTE ATLAS9 (IKuruca Il993l ) stellar 
atmosphere model with T eff = 24000 K, log(g) = 3.0 (cgs), turbulent velocity of 2 km s _1 
and solar metallicity. We fixed gravity brightening exponents and albedos to unity from 
theoretical values for stars at such temperatures. We assumed equal masses, circular orbits, 
and synchronous rotation, fitting for 7 parameters: the period P, time of primary eclipse 
r pr im(HJD), inclination i, T eff2 , the luminosity o f the primary in each band (L 1H ,Li K ) and 



the surface potential (Oj = Q2, see IWilsonlll979l . eq. 1). We defined convergence to be when 



the corrections for all adjusted parameters were smaller than their respective standard or 
statistical errors after three consecutive iterations. The best fit parameters are shown in 
Table [2j The ephemeris is 



prim 



2451775.102 ± 0.032 + 19.4513 ± 0.0011 x E (HJD). 



(1) 



A good fit required that the stars overfill their l i miting Roche lobes, which justi- 
fies using the overcontact mode of WD. [Martins et al.l (120061 ) adopted the largest filling 
factor allowed by NIGHTFA LL (1.3); we calculate a larger fill-out factor (as defined by 
Mochnacki fc Doughty! Il972l ) of F = 1.44. The critical surface potentials for the inner and 
outer surfaces under the above assumptions are fl in = 3.75 and Q out = 3.21. Assuming 
different values for the mass ratio q also produced good fits. However, a photometric mass 
ratio is not well constrained by our photometry; therefore , we did not a t tempt to solve for 
it. The inclination we derive is i — 71°, in agreement with iMartins et al.l (120061 ); however, if 
there is still unaccounted for blending, then the inclination angle could be larger. 

In Figures [2] and [3] w e show the H- and fC-band light curve model fit for IRS16SW. 
The IMartins et al.l (120061 ) data are plotted with unfilled symbols in Figure [3J Error bars 
for our data are set to 0.04 mag, corresponding to the typical variation seen for a constant 
star of similar magnitude (Peeples et al., in preparation). The fact that the data are fit 
well under the assumption of circular orbits is an indication that this assumption is sound; 
these data give no e vidence for an eccentric orbit. The eccentricity of e = 0.09 found by 
Martins et al.l (|2006h is derived from t heir radial velocity curve, which tends to yield nonzero 
eccentricities ( ILucy fc Sweeney! Il97ll ). Furthermore, the circularization time for a system 
with these physic al characteristices (discussed below) is only tens of thousands of years 
( 1Zahnlll975l . Il977l ). making it unlikely that we are observing IRS16SW pre-circularization. 



From the definition of Q, WD calculates the following best fit fractional radii for both 
stars in units of the orbital separation: the polar radius, r pole = 0.39; the radius in the 
plane of revolution and perpendicular to the line connecting the stars' center s, r^He = 0.41, 
and th e radius in the direction of L2, rb ac k = 0.47. The orbital separation of IMartins et al. 
hood ). (a x + a 2 )sinz = 132.8 ± 4.4i? , yields physical radii of R po i e = 54.5 ± 1.8 i? , 
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-Rside = 58.2 ± 1.9 R & , and -Rback = 62.7 ± 2.1 R Q . The WD visualization for this system is 
shown in Figure HI 



Discussion and Conclusion 



Using Kepler's law, iMartins et al.l (120061 ) find Mi ~ M 2 ~ 50M Q , placing the compo 
nents of IRS16SW among the most massive stars known. Until recently, the most massive 



stars measured in binaries were R136-38 in the Large Magellanic Cloud (5 7 My. lMassey et al. 



20021) and WR 22 (55 M , IRauw et al.lll996l . ISchweickhardt et all 11993 1. an evolved star in 
our Galaxy The current heavyweight cham pion is a Wolf-Rayet binary, WR 20a (82 & 83 
Mq, IRauw et al.l 12004 iBonanos et al.l 12004 ). in the young Galactic cluster Westerlund 2. 



The luminosity of IRS16SW poses a problem. Using a radius of R = 59.7R Q (the mean 
radius given by WD for an orb i tal se paration of 14O.6i? ) and an effective temperature of 
T e ff = 24400 K (jNajarro et al.lll997l ). we can estimate the luminosity L 



4nR 2 aT± s of 



each component as 4.4 x 10 erg s . The non-sphericity of IRS16SW will only drive this 
luminosity higher. For comparison, the Eddington luminosity of a 5OM star is L ec jd = 
1.3 x 10 38 (M/M Q ) = 6.5 x 10 39 erg s^ 1 . It is highly unlikely that each component of 
IRS16SW has been radiating s tably at nearly their Eddingt on luminos i ties f or eleven years 
(jliumphreys fc Davidson 1994 ) ; the combined photometry of Ott et al. Jl999h and this work 
span 1992-2002. Assuming the orbital separation as calculated by IMartins et al.l (120061 ) is 
correct (if it is smaller, then L/L p m will be even larger), this calculation implies that the 
temperature of lNajarro et al.l (119971 ) is an over-estimate. A change in the assumed T e g affects 
the WD model parameters; specifically, a decrease in T eS by a few thousand degrees Kelvin 
will decrease the inclination angle i, and thus increase the masses of the stars, by more than 
the formal la uncertainties given by WD. 

A radial velocity curve for the secondary is necessary to determine the value of the 
mass ratio q. It remains a puzzle as to why IRS16SW does not appear to be a double-line 
spectroscopic binary. It is readily apparent from the depths of the eclipses that the two 
stars have near- equal fluxes and fro m the depth ratio that they have near-equal surface 
brightnesses, yet IMartins et al.l (120061) see o n ly on e set of spectroscopic lines. However, since 
the spectral featu res used by IMartins et al.l (120061 ) are wind lines, with strong characteristic 
P Cygni profiles (jNajarro et al.lll997t l. differences in wind strength or small differences in 
the effective temperatures of the stars could easily conspire to make detection of the second 
set of lines difficult. 



We confirm that GCIRS 16SW is a massive eclipsing binary with both stars overflowing 
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their Roche lobes. We find a refined orbital period of 19.4513±0.0 011 days and an inclin ation 
of 71° with an assumed mass ratio of 1, supporting the findings by lMartins et al.l (120061 ) that 
the masses of the two stars are both ~ 5OM . The projected distance between IRS16SW 
and Sgr A* is 0.05 pc ~ 11000 AU (assuming a distance to the Galactic center of 7.6 kpc, 



Eisenhauer et al 



to Sgr A* flLu et al 



2005 ); in fact, IRS16SW is p art of a moving group that is likely bound 



20051 : IPaumard et al.l 120061 1 . As the lifetime of 50M© stars is ~ 4 Myr 



(ISchaller et al.l Il992l ). these observations are strong evidence that IRS16SW was formed 
within 0.1 pc of Sgr A* despite the tidal shear from the black hole which creates problems 
in star formation models. 
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Table 1. H- and K-band Photometry of GCIRS 16SW 



Band 


HJD 


Scaled 




-2450000. 


Magnitude 


H 


1769.6620 


9.630 


H 


1772.5868 


9.730 


H 


2542.5417 


10.041 


K 


2048.9194 


9.886 


K 


2051.7765 


9.527 



Note. - Both H- and fT-bands 
are scaled to have the same ampli- 
tude and mean magnitude as the 



Ott et all (119991) if -ban d data pre- 



sented by iMarfeisetaLl (J2006J). All 
errors are set to 0.04 mag, corre- 
sponding to the typical variation 
seen for a constant star of similar 
magnitude (Peeples et al., in prepa- 
ration) . Table CD is published in its 
entirety in the electronic edition of 
the Astrophysical Journal. A por- 
tion is shown here for guidance re- 
garding its form and content. 
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Table 2. Best-Fit Parameters From Combined H&.K Light-Curve Analysis 

With Wilson-Devinney Program 



Parameter Value 

Period, P 19.4513 ± 0.0011 days 
Time of primary eclipse, T prim 2451775.102 ± 0.032 

Inclination, i 70.85°±0.6° 

Temperature ratio, T 2 /Ti 0.96 

Surface potential, Q 3.51 

Light ratio in H, L 2 /Li 0.936 

Light ratio in K, L 2 /L x 0.939 

Radius, r pole) i = r po i ei2 0.39 

?"sidc,l = ^side,2 0.41 

?~back,l = ?"back,2 0.47 

Secondary temperature, T 2 23500 K 

Radius, -R po i e 54.5 ± 1.8 R Q 

-Rside 58.2 zh 1.9 Rq 

-Rback 62.7 ± 2.1 Rq 



Note. - First ten parameters are best-fit parame- 
ters from a combined HSzK light-curve analysis with 
the Wilson-Devinney (WD) program. The la uncertain- 
ties given by WD are unrealistically small, and thus not 
listed. The radii r pole , r side , and r back are in units of 
the orbital separation. The final four (physical) parame- 
ters are based on the orbital separation, (ai + a 2 ) sini = 
140.6 ± 4.7i? a nd the assumed effe ctive temperature of 
24400 K for T x ( TNajarro et al1ll997l ). 



9.6 9.8 10 10.2 

K magnitude 



Fig. 1. — H — K color residuals for DAOPhot photometry versus calibrated .fT-band magni- 
tude. No clear trend between color and magnitude is observed. The rms variation about a 
constant color is 0.05 magnitudes. The errorbar on the left shows the typical uncertainty in 
H — K color for a constant star of similar magnitude. See § [2] for further discussion. 
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Fig. 2. — Wilson-Devinney fit and residuals of an overcontact binary to the if-band light- 
curve of IRS16SW. The period is 19.4513 days; the model parameters have zero eccentricity 
and an inclination of 71°. The rms variation about the model fit is 0.06 magnitudes. 
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Fig. 3. — Wilson-Devinney fit and residuals of an overcontact binary to the i^-band light- 
curve of IRS1 6 SW. T he filled circles are the data presented here; the open squares are the 
Martins et al.l (120061 ) data. The period is 19.4513 days; the model parameters have zero 
eccentricity and an inclination of 71°. The rms variation about the model fit is 0.06 magni- 
tudes. 
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Fig. 4. — Wilson-Devinney visualization of IRS16SW at an orbital phase of 0.12. Axes are 
in units of R Q . 



